Hard diffractive quarkonium hadroproduction at high energies 
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We present a study of heavy quarkonium production in hard diffractive process by the Pomeron 
exchange for Tevatron and LHC energies. The numerical results are computed using recent exper- 
imental determination of the diffractive parton density functions in Pomeron and are corrected by 
unitarity corrections through gap survival probability factor. We give predictions for single as well 
as central diffractive ratios. These processes are sensitive to the gluon content of the Pomeron at 
small Bjorken-x and may be particularly useful in studying the small- a; physics. They may also be a 
good place to test the different available mechanisms for quarkonium production at hadron colliders. 
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I. INTRODUCTION 

Recently, it has been realized that much can be learnt 
about QCD from the wide variety of small- a; and hard 
diffractive processes, which are under intense study ex- 
perimentally. In particular, the diffractive heavy quarko- 
nium production has drawn attention because their large 
masses provide a natural scale to guarantee the applica- 
tion of perturbative QCD. There are several mechanisms 
proposed for the quarkonium production in hadron collid- 
ers [ll, II] , as the color singlet model, the color octet model 
and the color evaporation model. An important feature 
of these perturbative QCD models is that the cross sec- 
tion for quarkonium production is expressed in terms of 
the product of two gluon densities at large energies. This 
feature is transferred to the diffractive quarkonium pro- 
duction, which is now sensitive to the gluon content of 
the Pomeron at small- a; and may be particularly useful in 
studying the different mechanisms for quarkonium pro- 
duction. We notice that at hadron colliders the heavy 
quarkonium production is of special significance because 
they have extremely clean signature through their lep- 
tonic decay modes. 

Diffractive processes in hadron collisions are well de- 
scribed, with respect to the overall cross-sections, by 
Regge theory in terms of the exchange of a Pomeron 
with vacuum quantum numbers [d^- However, the na- 
ture of the Pomeron and its reaction mechanisms are 
not completely known. A good channel is the use of 
hard scattering to resolve the quark and gluon content 
in the Pomeron [J. Such a parton structure is natu- 
ral in a modern QCD approach to the strongly inter- 
acting Pomeron. The systematic observations of diffrac- 
tive deep inelastic scattering (DDIS) at HERA have in- 
creased the knowledge about the QCD Pomeron, provid- 
ing us with the diffractive distributions of singlet quarks 
and gluons in Pomeron as well as the diffractive struc- 
ture function 5]. In hadronic collisions, we shall charac- 
terize an event as single diffractive if one of the collid- 
ing hadrons emits a Pomeron that scatters off the other 
hadron. Hard diffractive events with a large momentum 
transfer are also characterized by the absence of hadronic 



energy in certain angular regions of the final state phase 
space (rapidity gaps). The events fulfilling the conditions 
of large rapidity gaps and a highly excited hadron rem- 
nant are named single diffractive in contrast to those in 
which both colliding hadrons remain intact as they each 
emit a Pomeron (central diffractive events). In the cen- 
tral diffractive processes, also known as double Pomeron 
exchange (DPE) processes, both incoming hadrons are 
quasi-elastically scattered and the final states system in 
the center region is produced by Pomeron-Pomeron inter- 
action. Here, we focus on the following single diffractive 
processes p + p{p) ^ p + J /"^ [T] -I- X and the central 
diffractive reactions, p + p{p) p + J/'^ [T] -I- p{p). 

Our starting point is the hard diffractive factorization, 
where the diffractive cross section is the convolution of 
diffractive parton distribution functions and the corre- 
sponding diffractive coefficient functions in similar way 
as the inclusive case. However, at high energies there 
are important contributions from unitarization effects to 
the single-Pomeron exchange cross section. These ab- 
sorptive or unitarity corrections cause the suppression of 
any large rapidity gap process, except elastic scattering. 
In the black disk limit the absorptive corrections may 
completely terminate those processes. This partially oc- 
curs in (anti)proton-proton collisions, where unitarity is 
nearly saturated at small impact parameters Q. The 
multi-Pomeron contributions depends, in general, on the 
particular hard process and it is called survival probabil- 
ity factor. At the Tevatron energy, ^/s = 1.8 TeV, the 
suppression (for single diffractive processes) is of order 
0.05-0.2 0, i, i, whereas for LHC energy, = 14 
TeV, the suppression appears to be 0.08-0.1 0, H, fiot - 
In the case of central diffraction, the theoretical esti- 
mates 0, Q give a survival probability factor between 
0.08 and 0.04 for Tevatron and LHC, respectively. There- 
fore, these corrections are quite important for the relia- 
bility of predictions for hard diffractive processes. 

The paper is organized as follows. In next section, we 
present the main formulas to compute the inclusive and 
diffractive cross sections for heavy quarkonium hadropro- 
duction. We also present the parameterization for the 
diffractive partons distribution in the Pomeron, extracted 
recently in DESY-HERA, and theoretical estimations for 
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the gap survival probability factor. In the last section we 
present the numerical results for Tevatron and perform 
predictions to future measurements at the LHC experi- 
ment. The compatibility with data is analyzed and the 
comparison with other approaches is considered. 



II. DIFFRACTIVE HADROPRODUCTION OF 
HEAVY QUARKONIUM 

Let us start by summarizing the formulas for inclu- 
sive and diffractive cross sections for quarkonium pro- 
duction in hadron colliders. For our purpose we will use 
the Color Evaporation Model (CEM) The main 

reasons for this choice are its simplicity and fast phe- 
nomenological implementation, which are the base for 
its relative success in describing high energy data (we 
quote Refs. P, Q for a review on CEM model and com- 
peting approaches). For completeness, we notice that 
the CEM model employed in the analysis it is subject 
of some criticisms (see Sec. 1.3.3 in Ref. [1]). In this 
model, the cross section for a process in which partons 
of two hadrons, hi and interact to produce a heavy 
quarkonium state, hi + h2 — > H{nJ^^) + X, is given 
by the cross section of open heavy-quark pair production 
that is summed over all spin and color states. All infor- 
mation on the non-perturbative transition of the QQ pair 
to the heavy quarkonium H of quantum numbers J^^ is 
contained in the factor F^jfc that a priori depends on 
all quantum numbers 

a{hi h2 H[nJ^P] X) = F„jPc (j{hi h2 ^ QQ X) , (1) 

where a{QQ) is the total hidden cross section of open 
heavy-quark production calculated by integrating over 
the QQ pair mass from 2niQ to 2mo, with mo is the 
mass of the associated open meson. The hidden cross 
section can be obtained from the usual expression for the 
total cross section to NLO. These hadronic cross sections 
in pp collisions can be written as 

app{y/s,mQ) = X! / dxidx2 ff{xi,^ip) f^{x2,fJ,p) 

X CTjj(Vs,'™q,A'f>mIj) , (2) 

where Xi and X2 are the fractional momenta carried by 
the colliding partons and f[ are the proton parton den- 
sities. The partonic cross sections are 
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with s the squared partonic center of mass energy, 
p — Am? / s and /j^'^'''' are the scaling functions given 
to NLO [l^. Here, we assume that the factorization 
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FIG. 1: (Color online) Total cross section {xp > 0) for J/\& 
hadroproduction (solid line) and charmonium states com- 
puted using CEM model. Single diffractive cross section, us- 
ing single- Pomeron exchange (dot-dashed line), and available 
accelerator data are also shown (see text). 



scale, fj,F^ and the renormalization scale, fif/, are equal, 
/i = = M-R- We also take /i = 2mQ, using the quark 
masses mc = 1.2 GeV and mi, = 4.75 GeV. These pa- 
rameters provide an adequate description of open heavy- 
flavour production [l^ . The invariant mass is integrated 
over 4m^ < s < 4m^ in the charmonium case and 
Ami ^ s < Am\ for T production. The factors F„jpc are 
experimentally determined [l3| to be Fn-- w 2.5 x 10~^ 
for J/* and Fn— « 4.6 x 10"^ for T. These coeffi- 
cients are obtained with NLO cross sections for heavy 
quark production p^ . 

In order to illustrate the success of the CEM model in 
describing current experimental data [Ij], the total cross 
section {xp > 0) for J/vP hadroproduction is numeri cally 
computed using MRST 2004 (NLO) set of partons [ij. 
The results are presented in Fig. [U where the curves 
are extrapolated up to the LHC energy, ^/s = 14 TcV. 
The agreement with the total cross section data is fairly 
good. For energies above Tevatron, possible corrections 
for gluon saturation may be sizable. The results can 
be directly compared to previous studies in charmonium 
production in Ref. [l^. The low x region is particularly 
relevant for J/^ production at the LHC as well as at 
Tevatron. For charmonium production, the gg process 
becomes dominant and information on the gluon distri- 
bution is of particular importance. The CEM model is 
able to address the common energy behavior of the dif- 
ferent quarkonium states. Following Ref. [15], in Fig. [T] 
the energy dependence of some charmonium states are 
also presented. We have used the experimental values 
a{'ip')/a{J/-^) ~ 0.14, cr(xc)/cr(J/V') 0.4 and the rel- 
ative weight a{'ip)/a{'tp') ~ 0.24. Here, we have denoted 
-0 by the directly produced IS cc state, in contrast with 
the experimentally observed J/4' (they can be originated 
also from radiative Xc decays). 



3 



In Fig. [21 the differential cross section ''■"p^^-"^ [^^^ 
is shown as a function of center of mass energy, -^/s. In 
this plot, the experimental data from RHIC (PHENIX 
point _16] at 200 GeV, from Au+Au electron measure- 
ments) and the higher energy data from Tevatron are 
also included. The physics involved in high energy ion- 
ion collisions might be very different from the one in 
proton-proton collisions for the case of PHENIX data 
and could explain the J/^ suppression due to a decon- 
fined medium. The Tevatron data indicates that correc- 
tion for gluon saturation may be important in bringing 
the theoretical curve closer to the experimental result. 
In particular, the typical values of Bjorken-x would be 
X ~ 2mcl\fs « 7 X 10"'' for Tevatron and x ~ 10"'' 
at the LHC taken at a relatively low momentum scale 
= 9 GeV^. In Fig. [3l the differential cross 



0' = = 

section B x "'"^ly''^ Iv^o is presented as a function of the 
center of mass energy (solid line). In addition, the mea- 
sured cross sections for the sum of the three T states 

a' _l_ _|_ x ") in the dilepton decay channel are shown 
. The agreement of the GEM model with accelerator 
data is very good and an extrapolation to the LHC is 
presented. The results give us considerable confidence 
in the extrapolation to the LHC energy. The predic- 
tions for direct bottomonium production in each state 
can be obtained by the ratios cr(T')/(T(T) ~ 0.36 and 
(T(T")/fT(T) ~ 0.27. 



A. Diffractive cross section - single— Pomeron 
exchange 

For the hard diffractive processes we will consider the 
Ingelman-Schlein (IS) picture Q, where the Pomeron 
structure (quark and gluon content) is probed. In the 
case of single diffraction, a Pomeron is emitted by one of 
the colliding hadrons. That hadron is detected, at least 
in principle, in the final state and the remaining hadron 
scatters off the emitted Pomeron. The diffractive cross 
section of a hadron-hadron collision is assumed to fac- 
torise into the total Pomeron-hadron cross section and 
the Pomeron fiux factor [4]. The single diffractive event, 
h\ + /i2 ^1,2 + i?[nJ'^^] -I- A, may then be written as 



dcf'^ {hi + hj + H[nJ^^] + X 
dx^d\ti\ 

Fnj^c X U/hX^^^^ M) 5- (IP + hj - 



gg + A),(3) 



where the Pomeron kinematical variable is defined 



as X ]p 



ergy in the Pomeron-hadron j system and ,Jsij = ^/s 
the center-of-mass energy in the hadron i-hadron j sys- 
tem. The momentum transfer in the hadron i vertex is 
denoted by i^. A similar factorization can also be applied 
to central diffraction, where both colliding hadrons can 
in principle be detected in the final state. The central 
quarkonium production, hi-\-h2^ hi + H[nJ'~^^\ + h2, is 



where 



is the center-of-mass en- 
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FIG. 2: (Color online) The differential cross section da/dy\y^o 
as a function of energy for the J/*!/ production (solid line). 
Single diffractive cross section (dot-dashed line) and acceler- 
ator data are also shown (see text). 



characterized by two quasi-elastic hadrons with rapidity 
gaps between them and the central heavy quarkonium 
products. The central diffractive cross section may then 
be written as. 



da^^ {hi + hj ^hi + H[nJ^^] + hj 



fwM4\\U\) U/A^w, M a{JP + lP^QQ + X). 

Here, we assume that one of the hadrons, say hadron 
hi, emits a Pomeron whose partons interact with par- 
tons of the hadron /i2. Thus the parton distribution 
xi,fi/hi{xi, fJ-'^) in Eq. ^ is replaced by the convolu- 
tion between a distribution of partons in the Pomeron, 
/3/a/i)(/3, /i^), and the "emission rate" of Pomerons by 
the hadron, /^/^(xp, i). The last quantity, fw/h{^w,t), 
is the Pomeron fiux factor and its explicit formulation is 
described in terms of Regge theory. Therefore, we can 
rewrite the parton distribution as 



dx^j^^ dx^^^ d\U\d\tj 



F, 



xifa/hi{xi, fJ, ) = J dxw J dp J dtfjp/h^xw, t) 

X /3/./p(/3,M')'5f/3-— V (4) 

Using the substitution given in Eq. ([1]), the hid- 
den heavy fiavour cross section can be obtained from 
Pomeron-hadron cross sections for single and central 
diffraction processes. 



d(7{w + h-.QQ + X) _ ^ fi/w{xi/x''^^ 



„(i) 



dxi dx2 — _ T- 
X fj/hAx2,HF)^ij{s,m'^Q,fi%,fijf)+ {1:=±2), (5) 
and 
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da (W + TP ^ QQ + X) 



dxi dx2 
fj/w {x2/x 



E 



fi/w (^xi/x^p^^ ; ^i^p 



In the numerical calculations, we will consider the 
difFractive pdf's recently obtained by the HI Collabora- 
tion at DESY-HERA [5]. The Pomeron structure func- 
tion has been modeled in terms of a light flavour singlet 
distribution consisting of u, d and s quarks and 

anti-quarks and a gluon distribution g{z). The Pomeron 
carries vacuum quantum numbers, thus it is assumed that 
the Pomeron quark and antiquark distributions are equal 
and flavour independent: q{p ~ q{p — 2j^^w, where 
is a Pomeron singlet quark distribution and Nf is the 
number of active flavours. Moreover, for the Pomeron 
flux factor, introduced in Eq. ([?]), we take the experi- 
mental analysis of the diffractive structure function [sf, 
where the x-jp dependence is parameterized using a flux 
factor motivated by Regge theory 0], 



fTP/-pixTp,t) = Aj 



2ap(t)-l 



(6) 



where the Pomeron trajectory is assumed to be linear, 
awit) = aff(O) -I- a'jpt, and the parameters Btp and 
are obtained from flts to HI FPS data IITI. 



B. Multiple-Pomeron exchange corrections 



where f2 is the opacity (or optical density) of the in- 
teraction. This quantity can be computed using a sim- 
ple one-channel eikonal model or more involved multiple- 
channel model. The opacity f2(s, 6) reaches a maximum 
in the centre of proton and becomes small in the periph- 
ery. Therefore, the survival factor depends on the spa- 
tial distribution of the constituents of the relevant sub- 
process. For instance, the spatial ^-distribution of sin- 
gle and double rapidity gap processes are assumed to be 
controlled by the slope B of the Pomeron-proton vertex, 
(3{t) (X exp(— _B|i|), and that there is no shrinkage coming 
from the Pomeron amplitude associated with the LRG 
in hard diffractive subprocesses. In addition, exp(— fi) is 
the probability that no inelastic soft interaction in the 
re-scattering eikonal chain results in inelasticity of the 
final state at energy s and impact parameter b. 

We will consider the theoretical estimates for < jS'p > 
from Ref. 10] (labeled KMR), which considers a two- 
channel eikonal model and rescattering effects. The sur- 
vival probability factor is computed for single, central 
and double diffractive processes at several energies, as- 
suming that the spatial distribution in impact parameter 
space is driven by the slope B of the pomeron-proton 
vertex. We will consider the results for single diffractive 
processes with 2B = 5.5 GeV~^ (slope of the electro- 
magnetic proton form factor) and without N* excitation, 
which is relevant to a forward proton spectrometer (FPS) 
measurement. Thus, we have < >|mr= 0.15, [0.09] 
and < >g^j^= 0.08, [0.04] for = 1.8 TeV (Teva- 
tron) [y/s = 14 TeV (LHC)]. There are similar theoretical 
estimates, as the GLM approach [19], which also consid- 
ers a multi-channel eikonal approach. We verify that the 
results are consistent with each other. We quote Ref. [l^ 
for a detailed comparison between the two approaches. 



In present calculation we will consider the suppression 
of the hard diffractive cross section by multiple-Pomeron 
scattering effects. This is taken into account through 
a gap survival probability factor, < \S\'^ >, which can 
be described in terms of screening or absorptive correc- 
tions [l^. This suppression factor of a hard process ac- 
companied by a rapidity gap depends not only on the 
probability of the initial state survive, but is sensitive 
to the spatial distribution of partons inside the incoming 
hadrons, and thus on the dynamics of the whole diffrac- 
tive part of the scattering matrix. The survival factor of 
a large rapidity gap (LRG) in a hadronic final state is 
the probability of a given LRG not be filled by debris, 
which originate from the soft re-scattering of the specta- 
tor partons and/or from the gluon radiation emitted by 
partons taking part in the hard interaction. Let A{s, b) 
be the amplitude of the particular diffractive process of 
interest, considered in the impact parameter, 6, space. 
Therefore, the probability that there is no extra inelastic 
interaction is 

^ I ^,2 ^_ Jd;'b\Ai.s,br e^p[-n{s,b)] 

J d^b\Ais,b)\^ ' 



III. RESULTS AND DISCUSSION 

Let us now present the predictions for hard diffractive 
production of heavy quarkonium. In the numerical cal- 
culations, the new HI parameterization for the diffrac- 
tive pdf's [1] has been used. The 'HI 2006 DPDF Fit 
A' is considered, with the cut xjp < 0.1. The single- 
Pomeron results are presented in Figs. [2] and [3] for J/^* 
and T production (dot-dashed curves), respectively. The 
single diffractive contribution is large, being of order 5- 
6 % from the inclusive cross section. The results cor- 
rected by multiple-Pomeron suppression factor are a fac- 
tor about 1/10 lower than the single-Pomeron ones. It 
should be stressed that sizable uncertainties are intro- 
duced by changing, for instance, quark masses and/or the 
renormalization scale. However, our purpose here is to 
estimate the diffractive ratios <t^ /atot, which are less sen- 
sitive to a particular choice. For sake of illustration, we 
have ^2^^|^£o = 85 (566) nb for the energy of = 2 
TeV and 159 (1770) nb for LHC energy, = 14 TeV. 
For T we have B ^^^|^=o = H (72) pb and 35 (386) pb 
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FIG. 3: (Color online) The differential cross section B x 
da/dy\y=o as a function of energy for the f _|_ f ' _|_ f " produc- 
tion (solid line). Single diffractive cross section (dot-dashed 
Une) and accelerator data are also shown (see text). 



for Tevatron and LHC, respectively. Numbers between 
parentheses correspond to values not corrected by sur- 
vival probability gap factor. 

Concerning the central difractive cross-sections, the 
predictions give small values. For instance, we have 
M^|CD^ = 18 nb and B'-^\% = 0.8 pb at = 2 
TeV. This gives little room to observe central diffrac- 
tive T events at the Tevatron but could be promising 
for the LHC. The values reach ^2^^|^i?o = 45 nb and 

d<y{r) iCD 



^^^^\v=o = 3 pb at = 14 TeV. Once again, the re- 
sults corrected by multiple-Pomeron suppression factor 
are a factor about 1/10 lower than the usual IS model. 
Let us now compute the diffractive ratios. We have de- 



fined the single diffractive ratio as i?sD 
the central diffractive ratio as i?cD 



.SD 



/da 
I dv 



dy 
/ dv 



and 
The 



dy 
da CD 
dy I dy^ 

results are summarized in Table I, where the diffractive 
ratios for heavy quarkonium production are presented for 
Tevatron and LHC energies. The multiple-Pomeron cor- 
rection factors are taken from KMR model. The numbers 
between parentheses represent the single-Pomeron calcu- 
lation. Based on these results we verify that the J/'^ 
and T production in single diffractive process could be 
observable in Tevatron and LHC, with a diffractive ratio 
of order of 1 % or less. This is a similar ratio measured in 
W and Z production at the Tevatron [23|. The predic- 
tions for central diffractive scattering are still not very 
promising, giving small ratios. However, the study of 
these events is worthwhile since their experimental sig- 
nals are quite clear. Finally, our theoretical prediction is 
in good agreement with the experimental measurement 



of CDF M, which found R, 



J/* 

SD 



1.45 ± 0.25 %. As the 



main theoretical uncertainty in determining the diffrac- 
tive ratio is the value for the survival factor, we consider 
the theoretical band 0.2lQ-0.15[l3| for Tevatron energy. 



TABLE I: Model predictions for single and central diffractive 
quarkonium production in Tevatron and the LHC. Numbers 
between parentheses represent the estimates using the single 
Pomeron exchange. 





Quarkonium 


7?SD (%) 


Rcu (%) 


2.0 TeV 
14 TeV 
2.0 TeV 
14 TeV 


J/* 

J/* 
(T + T' + T") 
(T + T' + T") 


0.93 (6.2) 
0.50 (5.9) 
0.78 (5.2) 
0.39 (4.3) 


0.2(2.5) 
0.15 (3.7) 
0.06 (0.7) 
0.03 (0.8) 



Therefore, our prediction gives i?iJ^eory ~ ^■^'^ ^ ^.19 %, 
which is consistent with the Tevatron determination. 

Our calculation can be compared to available litera- 
ture in diffractive heavy quarkonium production. For 
instance, in Ref. [l^l the large pt J/'^ production in 
hard diffractive process is computed using the color octet 
fragmentation mechanism and the normalized Pomeron 
flux [2^. They found a SD J/'^ cross section at large 
Pt (> 8 GeV) of order 10 pb and a diffractive ratio 
-RsD — 0.65 ± 0.15 %. This is closer to our calcula- 
tion for the SD ratio for Tevatron within the theoretical 
errors. Afterwards, in Ref. [13] the color-octet mech- 
anism combined with the two gluon exchange model(in 
LO approximation in QCD) for the diffractive J/^* pro- 
duction is considered. Now, the SD cross section is 
a{pp — > J/'^ X) — 66 nb. The comparison between 
these calculations shows the size of the large theoret- 
ical uncertainty. A related calculation, the J/^* + 7 
diffractive production, appeared in Ref. [25[ based on 
IS model (with normalized flux [235) ^'^'^ factorization 
formalism of NRQCD for quarkonia production. They 
found (j{pp [J/* + 7] X) = 3.0 pb (8.5 pb) and dif- 
factive ratio i?sD = 0.5 (0.15) % in central region at the 
Tevatron (LHC). These results are somewhat still com- 
patible with present calculation. 

Concerning central diffraction, there are some theoret- 
ical studies in literature. In Ref. [2^, the DPE pro- 
cess, p + p{p) ^ p + XJ + p{P)i is calculated using two- 
gluon exchange model in perturbative QCD. It is found 
the following DPE cross sections: cr(xco) — 735 nb and 
"■(Xfoo) — 0.88 nb. In the same work, it is found that 
'^"^^//+^^ = 2 nb/GeV and '^^f^ = 0.5 pb/GeV. Re- 

aydpT dyapT ' 

cently, in Ref. ;27'] the double-diffractive production of 
Xc and Xb mesons has been studied using also the Regge 
formalism and pQCD (including unitarity corrections). 
They found da/dy\y^o — 130 (340) nb for Xc production 
and da / dy\y=o — 0.2 (0.6) nb for Xb production for Teva- 
tron (LHC). It was verified that the exclusive produc- 
tion is a factor 10 larger than the inclusive rate. Finally, 
the diffractive x meson production bas been computed 
[2^ using the Bialas-Landshoff formalism and makes use 
of the DPEMC Monte-Carlo simulation for small-mass 
diffractive production. They found iicD — 6. 5 (1.6) % 
and R^\^ = 22 (1.83) % for Tevatron (LHC). In summary, 
the theoretical predictions for exclusive meson produc- 
tion are still quite distinct and more detailed studies are 



6 



deserved. 

Concerning the central difFractive quarkonia produc- 
tion, there are some competing processes that have been 
studied recently in literature. The first one is the ex- 
clusive pp p + M + p process, with M — J/^* (T), 
mediated by photon exchange. In its simplest version, 
this process can be computed using Weizsacker- Williams 
approximation where the proton is treated as a source 
of quasi-real photons that collide with the other pro- 
ton and produce vector mesons. For instance, it has 
been treated in detail (going beyond the equivalent pho- 
ton approach) in Refs. [H, [sO, HH. Here, I quote the 
central scenario presented in recent study of Ref. [30| : 
da-y/dy {pp ppJ/'^)\y=o ~ 5 (15) nb and da-y/dy {pp 
ppT)\y^Q ~ 5(31) pb for Tevatron (LHC) energy. Sim- 
ilar results are obtained in Ref. [31|. Another pro- 
cess leading to exclusive meson hadroproduction is the 
Pomeron-Oddcron fusion mechanism, which can be in- 
vestigated in the framework of Rcggc theory (see, for ex- 
ample Ref. [S^]) or in the framework of fc^-factorization 
[30| (the mass of the heavy vector meson supplies the 
hard scale in the process, allowing a treatment within 
perturbative QCD). For the QCD Odderon we quote re- 
sults in REf. 30], where daodd/dy {pp — > ppJ/'i')\y=o — 
1.3(0.9) nb and daodd/dy {pp ppT)\y=o — 4(5) pb 
for Tevatron (LHC) energy. Let us compare these val- 
ues with the results for central Pomcron-Pomeron pro- 
duction presented here: dacD/dy {J/'if)\y=o — 18(45) 
nb and dacD / dy {T)\y=o ~ 32 (120) pb for Tevatron 
(LHC) energy, which include also the absobtive correc- 
tions. This indicates that the central diffractive contribu- 
tion is substantially larger than photon and/or Odderon 
exchange. However, the different competing processes 
could be identified by using additional experimental cuts 
as the outgoing proton momenta distibution. 

Finally, some comments are in order on the limita- 
tions of the present study. The theoretical uncertain- 
ties in both, quarkonium production and diffraction, are 
large and changes in the quarkonium models are likely to 
be compensated by changes in the diffractive part, the 
value of the heavy quark mass and of the factorization 
scale (and vice versa). A more detailed study should be 
done to clarify these issues. For example, in order to 
demonstrate the sensitivity to the diffractive gluon, it 
would have been nice to propagate the uncertainty band 
of the HI Fit A gluons or to compare with the Fit B 
gluons. Also the variation of the scales and heavy quark 



mass would have been easily possible and interesting for 
the absolute cross sections in Figs. 1-3 but also for the 
diffractive ratios to see how much these variations con- 
tribute to the theoretical error which is dominated by the 
uncertainty of the gap survival factor. Reliable estimates 
of the uncertainties are of particular importance because 
the existing theoretical predictions for exclusive meson 
production are quite distinct. 

In summary, we have presented predictions for diffrac- 
tive heavy quarkonium production at the Tevatron and 
the LHC. We use Regge factorization, corrected by uni- 
tarity corrections modeled by a gap survival probabil- 
ity factor (correction for multiple-Pomeron exchange). 
The perturbative formalism for meson hadroproduction 
is based on the CEM model, which is quite successful in 
describing experimental results for inclusive production. 
For the Pomeron structure function, recent HI diffractive 
parton density functions have been used. The results are 
directly dependent on the quark singlet and gluon con- 
tent of the Pomeron. We estimate the multiple interac- 
tion corrections taking the theoretical prediction a multi- 
channel model (KMR) , where the gap factor decreases on 
energy. That is, < |5p >~ 15 % for Tevatron energies 
going down to < 15*^ >~ 9 % at LHC energy (for single 
diffractive process) . We found that at the Tevatron single 
and central diffractive J/^l/ and T production is observ- 
able with a single diffractive ratio i?^'-' (Tevatron) that is 
between 1 % (charmonium) and 0.5 % (bottomonium), 
with lower values at the LHC. In particular, we predict 
-^fhoory ~ 1.12±0.19%, which is in agreement with Teva- 
tron measurements. The central diffractive cross sections 
for quarkonium production are still feasible to be mea- 
sured, despite the very small diffractive ratios. In this 
case, the theoretical model dependence is still very large 
and detailed studies are deserved. 
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